The charmed and bottom meson spectrum from lattice NRQCD 



Randy Lewis a and R. M. Woloshyn b 

a Department of Physics, University of Regina, Regina, SK, Canada S4S 0A2 
b TRIUMF, 4004 Wesbrook Mall, Vancouver, BC, Canada V6T 2A3 

The mass spectrum of S and P-wave mesons containing a single heavy quark has been computed using quenched 
lattice nonrelativistic QCD. Numerical results have been obtained at first, second and third order in the heavy 
quark expansion, so convergence can be discussed. The computed spectrum of charmed and bottom mesons is 
compared to existing model calculations and experimental data. 



1. MOTIVATION 



lations are compared to theoretical models. 



The masses of the ground state doublets of S- 
wave heavy-light mesons (D^*\ Dg, and 
Bs ) are well known experimentally, and thus 
provide a testing ground for the lattice method. 
In particular, they can be used to examine the 
convergence properties of the heavy quark expan- 
sion in nonrelativistic QCD (NRQCD), Q] which 
is the most economical lattice action available for 
heavy quark physics. Q 

The spectrum of P-wave heavy-light mesons is 
not yet understood, experimentally nor theoreti- 
cally, so here the lattice can offer predictions. Do 
the masses follow the ordering that is observed 
in the hydrogen atom? What are the magnitudes 
of the mass splittings among P-wavcs tens of 
MeV or hundreds? 

Lattice studies of the P-wave spectrum have 
been performed for an infinitely- heavy quark Q. 
A subset of the P-wave states have also been ob- 
served by using a relativistic (clover) action in 
the charmed region, and extrapolating to the bot- 
tom region Recently, the lattice NRQCD ap- 
proach has been used to work directly with B** 
mesons as well as D** mesons [||J7|. 

In the following sections, the lattice NRQCD 
method is briefly reviewed and the results of nu- 
merical simulations are presented. Implications 
for the convergence of the lattice NRQCD expan- 
sion are discussed, and predictions for the B** 
and D** mass spectra from various lattice simu- 



2. LATTICE NRQCD METHOD 

The application of lattice NRQCD to heavy- 
light mesons entails an expansion of the action in 
inverse powers of the heavy quark mass M , 

S = S q . G + J d 4 x Qt (iD t + ^^SH^j Q, (1) 

where S* 9 ,g has no heavy quark fields, and||,^| 
SH = SH^ + 5H {2) + 5HW + 0(1/Af 4 ), (2) 
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The chromoelectric and chromomagnetic fields 
are Ei — Fa and Bi — eijkFijk/2 respectively, 
and A denotes the lattice derivative. The coef- 
ficients Ci are usually set to their classical values 
(c, = 1) in lattice NRQCD simulations. This as- 
sumes that perturbative corrections at momen- 
tum scales beyond the lattice cutoff (i.e. the in- 
verse lattice spacing) are small. Nonperturbative 
contributions to the a coefficients are not ne- 
glected, and appear explicitly in Eqs. (|f|-||) via 
the tadpole factors, mentioned below. 

A tilde on any quantity indicates that the lead- 
ing classical discretization errors have been re- 
moved. Besides this classical improvement, a 
mean-field form of quantum improvement is also 
implemented through the introduction of two 
"tadpole factors", u s and u t , following the well- 
known suggestion of Lepage and Mackenzie ||. 
Notice that the action of Eqs. (jf]||) allows dif- 
ferent lattice spacings, a s and a t , in the spatial 
and temporal directions; hence two tadpole fac- 
tors are required. 

Various choices for the light quark and gauge 
field action, S q ,G of Eq. (jl|), have been considered. 
In the lattice NRQCD simulations of Refs. [QjH, 
the gauge field portion is classically and tadpole- 
improved, whereas the authors of Refs. HQ chose 
the unimproved Wilson gauge field action. For 
the fcrmion portion, Refs. pjM fl] employ the 
tadpole-improved clover action|K ] while Rcf . |5|| 
uses a more highly improved D234 action [[nf. 
When the heavy quark mass is sufficiently large 
and the lattice spacing sufficiently small (but not 
small compared to the inverse heavy quark mass), 
one expects the predictions of these actions to 
agree. 

3. S-WAVE HEAVY-LIGHT MASSES 

Lattice NRQCD simulations typically produce 
a clear signal for the S-wave heavy-light meson 
masses, so they are convenient observables to use 
for testing the heavy quark expansion. Since 
the leading contribution to the heavy quark mass 
does not appear explicitly in the NRQCD action, 
only mass differences (or more generally energy 
differences) can be computed. The absolute mass 
scale can be extracted from the kinetic energy of 
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Figure I. S-wave splittings of charmed mesons 
at first, second and third order in NRQCD from 
Ref. |5). Solid lines simply connect related data 
points; they are not fits. 
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to the desired order in l/M^in, where E p and 
Eq are taken from lattice simulations, and the 
"kinetic" mass Mki n is thus obtained. 

The first three orders in the 1 /M expansion of 
Eqs. (|f|-^) have been studied in Refs. p]J| for two 
different light quark actions. The charmed results 
of Ref. || are displayed in Fig. [l| In this plot, 
subleading corrections to the D s — D splitting are 
insignificant, and the physical result agrees well 
with the experimental result pjfl, 
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The D*-D splitting of Fig. | has 0{l/M 2 ) and 
0(1/ M 3 ) corrections which are roughly equal, at 
about 20% each, and the D* — D s splitting ac- 
quires 20% and 10% corrections at 0(1/M 2 ) and 
0(1 /M 3 ) respectively. Firm statements about 
convergence, or lack thereof, for the 1/M ex- 
pansion are difficult to support with only these 
data. The first correction terms are not far 
from their expected size, 0(Aqcd/M charm), but 
the 0(1 /M 3 ) corrections are perhaps larger than 



Table 1 

Quenched lattice NRQCD results for the S- 
wave bottom mesons Q compared to experimen- 
tal data|J|], in units of MeV. 

quenched experiment 



B* 



B 
B 
B 



92±3 
25±2 
26±1 



90±2 

46 
47±3 



might have been hoped. These spin splittings are 
clearly smaller than the experimental values fll2|| 
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as typical of quenched lattice calculations. |ljj4j- 

The 1 /M expansion is nicely suited to the bot- 
tom mesons and, even with 2000 gauge field con- 
figurations, the simulations of Ref. [|| found the 
nonleading contributions to be smaller than sta- 
tistical uncertainties. The quenched results are 
shown beside the experimental data in Table [l]. 

4. P-WAVE HEAVY-LIGHT MASSES 

It is more difficult to extract a signal for the P- 
wave masses from lattice simulations. Creation 
operators which extend over more than a single 
lattice site need to be tuned to enhance the signal. 

In Ref. §, the P-wave masses were studied 
at 0(l/M), 0(1/M 2 ) and 0(1/M 3 ) for both 
charmed and bottom mesons. Subleading contri- 
butions were found to be smaller than the statis- 
tical uncertainties. Fig. |^ compares the splitting 
between two P-wave masses (J=0 and J=2) as ob- 
tained from lattice calculations and from various 
models. 1 13 For the B** system, most methods 



predict a splitting between zero and 80 MeV. 

The discrepancy between the two lattice re- 
sults is somewhat disconcerting. Both compu- 
tations employed quenched lattice NRQCD and 
both had a temporal lattice spacing of about 
0.1 fm. Ref. ^] used an improved gauge action 
with a spatial lattice spacing of 0.2 fm, a highly- 
improved (D234) light-quark action and an en- 
semble of 2000 configurations. Ref. Q used the 
unimproved (Wilson) gauge action with a spatial 
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(b) B 2 * - B Q * and B s2 * 
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Figure 2. Splittings between the Pq and Pi 
masses for the D** and B** systems. Open(solid) 
symbols involve an s(u,d) quark. 



lattice spacing near 0.1 fm, the improved clover 
action for light quarks and an ensemble of 102 
configurations. 

In Ref. §], the ratio of P 2 /P correlation func- 
tions falls monotonically as Euclidean time in- 
creases, and of course statistical uncertainties 
grow as Euclidean time increases. Thus, if the 
meson mass is extracted from a region which be- 
gins too near the initial time, the predicted mass 
splitting would be too large. With this in mind, it 
is noted in Ref. || that, even at the smallest Eu- 
clidean times, those lattice data cannot reproduce 
the large splitting of Ref. |6). Instead, Ref. || 
arrives at a plateau-independent upper bound of 
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Figure 3. S-P splittings for heavy-light mesons. 
Open(solid) symbols involve an s(u,d) quark. 



100 MeV for the B^-B^ splitting. Lattice efforts 
are ongoing, and we look forward to a clearer un- 
derstanding of the D** and B** spectra. 

In Fig. |H the splitting between S and P-wave 
mesons is shown. The lattice results of Boyle are 
interesting in that they do not rely on NRQCD. Q 

5. SUMMARY 

The convergence of lattice NRQCD has been 
studied up to 0(1/M 3 ) for both S and F-wave 
heavy-light mesons. The only cause for concern 
seems to be the S-wave spin splitting for charmed 
mesons, where a definitive statement cannot yet 
be made about convergence. 

The P-wave masses have recently come under 



study with lattice NRQCD. The results of Ref. g 
indicate that the D** splittings are O(50 MeV) or 
smaller, while the B** splittings are O(10 MeV) 
or smaller. However, the lattice work reported 
in Ref. Q is not consistent with this, so further 
studies will be instructive. 
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